Circulating Nucleic Acid-Based Biomarkers of Type 2 Diabetes

Type 2 diabetes (T2D) is a deficiency in how the body regulates glucose. Uncontrolled T2D will result in chronic high blood sugar levels, eventually resulting in T2D complications. These complications, such as kidney, eye, and nerve damage, are even harder to treat. Identifying individuals at high risk of developing T2D and its complications is essential for early prevention and treatment. Numerous studies have been done to identify biomarkers for T2D diagnosis and prognosis. This review focuses on recent T2D biomarker studies based on circulating nucleic acids using different omics technologies: genomics, transcriptomics, and epigenomics. Omics studies have profiled biomarker candidates from blood, urine, and other non-invasive samples. Despite methodological differences, several candidate biomarkers were reported for the risk and diagnosis of T2D, the prognosis of T2D complications, and pharmacodynamics of T2D treatments. Future studies should be done to validate the findings in larger samples and blood-based biomarkers in non-invasive samples to support the realization of precision medicine for T2D.


Introduction
Type 2 diabetes (T2D) remains a significant clinical burden worldwide. Approximately 462 million individuals (6.28% of the world's population) lived with T2D in 2017, and the prevalence continues to rise [1]. Better screening, diagnosis, and treatment approaches are needed to combat T2D. Early screening would ensure the timely implementation of lifestyle interventions for those at risk of developing T2D [2]. Early diagnosis is substantially beneficial for T2D patients. Even though T2D incidence peaks at age 55-59, many develop the disease at earlier ages and are often undiagnosed for years [1,3]. Early diagnosis would allow early treatment and prevent T2D complications, which are more challenging to treat [4]. Personalized T2D management, where drug choice is based on each patient's characteristics, would maximize glycemic control efficiency and minimize side effects [5].
Numerous studies have been conducted to identify T2D biomarkers. Several reviews have summarized genetic and non-genetic susceptibility/risk biomarkers for developing T2D [15][16][17]. Another review summarized omics-based biomarkers for diagnosing different types of diabetes and included non-invasive biomarkers from urine and saliva [18]. Research on omics-based and non-invasive T2D biomarkers has increased in recent years and proposed potential novel biomarkers. So, this review aims to summarize candidate human T2D biomarkers from a genomic, transcriptomic, and epigenomic perspective. We looked into T2D studies reporting all biomarkers categories, including risk and diagnostic biomarkers from blood, urine, and saliva samples. The gene names and symbols are listed in Table S1 (see Supplementary Materials). Predictive A biomarker used to identify individuals who are more likely than similar individuals without the biomarker to experience a favorable or unfavorable effect from exposure to a medical product or an environmental agent BRCA1/2 mutations can identify ovarian cancer patients likely to respond to PARP inhibitors Pharmacodynamic/response A biomarker used to show that a biological response has occurred in an individual who has been exposed to a medical product or an environmental agent HbA1c may be used to assess diabetes control after treatment Safety A biomarker measured before or after an exposure to a medical product or an environmental agent to indicate the likelihood, presence, or extent of toxicity as an adverse effect Neutrophil count can be used to adjust dose for patients on cytotoxic chemotherapy.

Genomics Biomarkers for T2D
The science of genomes, or "genomics", studies the structure, expression, and function of the whole DNA sequences of an organism. It has rapidly expanded towards a more functional level, studying the evolution, mapping, and editing of genomes [19]. Genomics analyzes genetic variants, such as single nucleotide polymorphisms (SNPs) and chromosomal abnormalities related to medical conditions [18]. T2D genomics studies used TaqMan qPCR assays, next-generation sequencing, and DNA microarray technologies.

Blood
Among all sample types, blood biomarkers represent the most accessible and most studied, given their properties and ease of collection in clinical practice [20]. A wide variety of T2D-related biomarkers in blood using genomic data has been reported ( Table 2).
Heritability is a strong predictor of T2D (26-69% depending on the age of onset), thus motivating the search for genetic predictors of T2D [21]. One way to translate genetic data to a predictive measure of disease susceptibility is to add the risk effects of loci into a polygenic risk score (PRS) [22,23]. Prediction accuracy of a PRS is often assessed by measuring the area under the receiver operating characteristic curve (AUC). The AUC compares the rates of true positives and false positives, accounting for the overall performance of predictive models [24]. The first PRS for T2D was developed in 2006 using three genetic variants: KCNJ11, PPARG, and TCF7L2 [25]. Two years later, three PRSs were developed using more SNPs. Two of them had 18 SNPs, and the third one had 11. Their AUC was 78%, 90%, and 74%, respectively [26][27][28]. In 2014, a 62-SNP PRS was developed using the DIAGRAMv3 Panel of Genes [29], with an AUC of 90% [30]. In 2016, Chikowore analyzed a South African population; a PRS using only four SNPs was created, taking into account sex, age, BMI, and systolic blood pressure as clinical risk factors, with an AUC of 66.5%. This research was the first to develop PRS in non-caucasian ethnic group with a high sensitivity and specificity [31]. Recent PRS models include many more SNPs, thanks to the larger genome-wide association studies (GWAS) in recent years. Furthermore, in 2016, Lall and colleagues published a PRS using 1000 SNPs from the DIAGRAM Panel of Genes, with an AUC of 77% [32].
Besides using SNPs in a PRS, some studies demonstrate other ways of utilizing SNPs as risk biomarkers. Ding and colleagues demonstrate that a low circulating level of sex hormone-binding globulin is a potential risk biomarker of T2D. Nevertheless, the clinical use of SHBG needs further examination [33]. In 2015, a study showed that endogenous bilirubin and the associated SNP, rs6742078 in the UGT1A1 locus, are associated with the risk of T2D, making it a candidate as a risk biomarker [34]. Two years later, Wheeler identified 60 common genetic variants associated with HbA1c using genome-wide association metaanalyses from 82 European, African, East Asian, and South Asian ancestry cohorts. They found out that a G6PD deficiency can be clinically silent until illness strikes. Therefore, screening with direct glucose measurements in people with G6PD deficiency may be helpful as a risk predictor of T2D [35].
One of the first studies researching predictor biomarkers of T2D was done in 2010. In the study, the levels of plasma GAPDH, representing total cell-free DNA, were measured. The levels of cell-free GAPDH were significantly higher in the plasma samples of T2D diabetic patients, becoming a biomarker candidate [36]. The most recent study was done in 2019, where they assessed the potential role of mitochondria in T2D by analyzing bloodbased indices of mitochondrial DNA copy number (mtDNACN) and cell-free mitochondrial DNA (CFmtDNA). As a result, they found a significant locus in the LRRK2/MUC19 region, indicating that mitochondrial dysfunction is intimately linked to T2D prediction, therefore being a candidate biomarker [37].

Urine
Urine samples have several advantages above blood: the collection is easy and noninvasive, and the samples are available in large volumes. Anyone can collect urine samples, unlike blood samples which require clinical personnel [38]. Urine is not generally used for DNA analysis because the extracted DNA is less stable than blood. So, most T2D genomics studies collect blood and urine samples to study the associations between blood genetic variants and urine proteins (Table 3). However, recent studies have improved DNA extraction from urine [39]. One recent T2D study reported associations between mitochondrial DNA (mtDNA) from urine samples with T2D, demonstrating the possibility of urine-extracted DNA as biomarkers for T2D [40].  Blood (plasma or serum) and urine qRT-PCR 4668, 0, 2290 ** Urine albumin-to-creatinine ratio and WFS1 (rs10010131) [41] Blood (plasma) and urine qRT-PCR 0, 290, 285 ** Urine albumin and PGC-1α (rs8192678) [42] Blood and urine qRT-PCR 35, 0, 42 ** Urine creatinine ratio with mtDNA [40] ** T2D patients with complications or comorbidities.
All T2D genomics urine studies looked into diabetic kidney disease (DKD), including diabetic nephropathy (DN) ( Table 3). DKD is one of the significant complications of T2D and is the leading cause of end-stage renal failure. About one-third of diabetes patients will develop kidney disease, so identifying DKD biomarkers is essential for early prevention [43]. T2D studies on urinary biomarkers usually include patients with normo-, micro-, and macroalbuminuria.
Two studies found associations between SNPs measured from blood samples with urine proteins. A prospective study reported that rs10010131 in the WFS1 gene is associated with a higher estimated glomerular filtration rate (eGFR) in T2D patients with increased albuminuria. This result suggests a potential role of WFS1 as a risk biomarker of T2D and its kidney complications [41]. The second study examined the association between urine albumin and rs8192678 in the PGC-1α gene, a genotype previously related to nephropathy in T2D patients. The genotype was found to be associated with a 70% increase in urine albumin excretion in T2D patients with proteinuria compared to T2D patients with normoalbuminuria, making PGC-1α a candidate prognostic biomarker [42].
The only T2D study on urine-extracted DNA looked into mtDNA from healthy controls and T2D patients with and without proteinuria. Plasma and urine mtDNA content significantly differed between T2D patients and controls, where a reduction in plasma mtDNA content and increased urinary mtDNA/creatinine ratio were observed. The study did not report significant differences between T2D patients with and without proteinuria. So, the findings suggest that mtDNA could be a diagnostic biomarker of T2D [40].

Other Non-Invasive Biomarkers and the Use of Metagenomics
Saliva is another non-invasive sample to research biomarkers using genomic technology. Saliva-extracted DNA has comparable quantity and quality with blood, making saliva a suitable source of DNA for genetic studies. So, unlike urine, saliva has been used in large genetic epidemiological and metabolic studies [44]. However, there are no T2D studies that profile human DNA from saliva samples. One study profiled polymorphisms for CHGA from peripheral blood leukocytes and measured CHGA concentration in saliva ( Table 4). The study found associations between two polymorphisms (rs9658635 in the promoter region and Glu264Asp in exon 6) with higher salivary CHGA production [45]. It remains to be seen whether these DNA biomarkers can also be detected in saliva samples.
In addition to human DNA, microbial DNA can also be extracted from saliva. Metagenomics studies the DNA sequences of an entire community of microorganisms [46]. The metagenomics studies aim to find associations between microbial profiles and human phenotypes. The technologies most frequently used in this area are 16S rRNA sequencing and whole-genome sequencing [18]. A recent study reported a higher relative abundance of Bulleidia, Ruminococcaceae, and Helicobacter pylori in T2D patients compared to healthy controls. However, the study was done in only nine T2D patients, so future studies are needed to confirm the findings [47] (Table 4). T2D genetic studies have also been done using fecal samples. The genetic material from fecal samples is mostly bacterial, although human DNA can be detected in small amounts. There are efforts to improve human DNA extraction from stool samples, but significant challenges still exist for its use in population studies [53]. Therefore, all T2D studies with fecal samples are metagenomics studies (Table 4).
In 2016, a study showed evidence that the inflammation of the gut increased the values of biomarkers related to T2D. Microbial signatures of Akkermansia muciniphila demonstrate the existence of inflammation, increasing the risk of T2D before they are reflected by clinical markers [49]. Other metagenomics studies identified Akkermansia muciniphila, Ruminococcus torques, Ruminococcaceae, Lachnospiraceae, and Enterobacteriaceae to be associated with T2D [54][55][56][57]. In a critically high-risk population, a recent study demonstrated associations between Escherichia, Veillonella, Blautia and Anaerostipes with T2D [51]. Another study looked into the abundance of antibiotic resistance genes (ARG) in fecal microbiome profiles. The study found ARG enrichment in T2D patients, which could be potential biomarkers of T2D [50]. It remains to be seen if these T2D-associated changes in the gut contribute to T2D pathogenesis.

Transcriptomics Biomarkers of T2D
Transcriptomics is the comprehensive study of total gene expression levels in a cell or organism [58,59]. Most T2D messenger RNA (mRNA) research does not qualify this definition, as most studies profiled specific genes related to the study population. Nevertheless, the studies have identified potential mRNA biomarkers of T2D. The most common profiling method for T2D mRNA studies is targeted quantitative reverse-transcription PCR, where isolated RNAs are reverse transcribed to cDNA, then the expression levels of specific genes are measured quantitatively using qPCR. A few more recent T2D studies use microarray and RNA sequencing to measure many genes, but only in a small number of samples. Then, the top differentially expressed genes are validated using qRT-PCR in a larger number of samples. Currently, most T2D mRNA studies have been conducted with blood, but there is an increasing number of studies on non-invasive samples, such as urine and saliva.

Blood
Blood mRNA T2D studies have profiled gene expression in many sample types, from whole blood samples to leukocytes (Table 5). Only one study compared gene expression levels between two sample types: serum and urine [60]. There is little consensus on mRNA biomarkers because all studies selected genes related to their study population.
Blood mRNA studies on T2D patients without complications have identified elevated levels of inflammatory and senescence markers, as well as metabolically relevant receptor genes [61][62][63][64]. Plasma NF-kB expression levels were increased in pre-diabetic patients and further increased in T2D patients [62], suggesting that NF-kB can be a biomarker of early T2D diagnosis. In contrast, blood leptin mRNA levels were reduced in T2D patients compared to non-diabetic controls [65]. The discriminatory power of leptin mRNA was comparable to fasting blood glucose and HbA1c in distinguishing T2D patients and non-diabetic controls (AUROC = 0.95, 95% CI: 0.89-0.98) [65].
Studies on T2D patients with complications have focused on profiling genes related to T2D, inflammation, or specific diabetic complications [60,[66][67][68][69][70][71]. For example, a study on diabetic kidney disease (DKD) patients selected AEBP1 from a set of estimated glomerular filtration rate (eGFR)-correlated genes [66]. Another study on T2D patients with bone mass reduction or osteoporosis profiled BSP mRNA, a marker of bone metabolism [71]. The overall findings are increased inflammation and reduced expression of protective genes in patients with T2D complications.
One study identified a circulating biomarker for treatment response. A clinical trial administered placebo or zinc supplementation in overweight T2D patients and profiled SOD1 expression in leukocytes [72]. The mRNA levels were significantly elevated in patients who received supplementation, compared to placebo and levels before treatment, and correlated with serum zinc concentrations. So, blood mRNAs can also be response biomarkers to evaluate treatment in T2D patients. More studies are needed to validate and identify other response biomarkers for T2D patients.

Urine
Similar to genetics studies, urinary mRNA studies have compared T2D patients with and without DKD ( Table 6). The first study looked into the expression levels of TTP and inflammatory markers IL-6 and IL-18. TTP encodes tristetraprolin, which is an anti-inflammatory protein. A significantly reduced TTP expression level was found in T2D patients, and the levels are further reduced in patients with macroalbuminuria. In contrast, IL-6 and IL-18 were significantly elevated in T2D patients, especially those with complications [60].
Another study looked into podocyte-associated markers: NPHS1, NPHS2, and PODXL [73]. Renal podocyte injury is pathologically significant in DN progression, so podocyte-associated genes might serve as monitoring biomarkers to evaluate kidney health. Urinary podocyte marker levels were significantly elevated in all T2D patients, including those with complications. Patients with macroalbuminuria had higher podocyte marker levels than normoand microalbuminuria patients. These markers distinguished microalbuminuria patients from those without complications in univariate analyses, with AUCs of 0.685-0.961 [73].
A study with 242 participants used the presence of albuminuria and eGFR to define their study groups [74]. This study evaluated gene expression in urinary exosomes using RNA sequencing and validated eight mRNAs (UMOD, SLC12A1, NDUFB2, OAZ1, PPARGC1A, NFE2L2, CD24, SMAD1) using qRT-PCR. So, this is the first agnostic urinary mRNA study where markers were not selected from previously available knowledge. Among the eight mRNAs, UMOD and SLC12A1 levels were elevated in T2D patients with complications compared to those without complications and healthy controls. Significant increases were also seen in NDUFB2 and OAZ1 in T2D patients without complications compared to healthy controls. The four mRNAs were significantly correlated to HbA1c%, but only OAZ1 was correlated to T2D duration. A multi-gene signature was also developed using the eight mRNAs and linear discriminant analysis with cross-validation. The classifier distinguished severe DKD patients from T2D with an AUC of 0.90 [74].

Other Non-Invasive Biomarkers
Unlike genetic T2D biomarkers, there is limited research for mRNA biomarkers from saliva samples. One study reported differential salivary gene expression from 13 T2D patients: elevated KRAS, SAT1, SLC13A2, and TMEM72 and reduced EGFR and PSMB2 gene expressions [75]. The six biomarkers were first identified through microarray, then validated in 13 T2D patients and 13 healthy controls. A logistic model with four out of the six salivary biomarkers (KRAS, SAT1, EGRF, and PSMB2) discriminated T2D patients from healthy controls (AUC: 0.917, 95% CI: 0.809-1.000). So, salivary mRNAs are promising T2D biomarkers, but more studies on these samples are needed.

Epigenomics Biomarkers for T2D
Epigenomics studies the whole gene regulation of an organism or cell, including DNA methylation and microRNAs (miRNAs) [59,76]. Epigenetic factors regulate gene expression without modifying DNA sequence. For example, DNA methylation of CpG islands influences transcription factor binding, resulting in stable silencing of gene transcription [77]. Gene expression can also be regulated post-transcriptionally through miRNAs, which bind to target mRNAs and typically lead to translational repression [78]. The miRNA-mRNA interactions can be binary off-switch interactions or tuning interactions in response to an environmental cue [79].
There is growing interest in identifying circulating epigenetic biomarkers for T2D, especially since T2D has a substantial environmental influence, such as sedentary lifestyle and obesity. Identifying risk, predictive, monitoring, and response biomarkers are among the most common research aims in T2D epigenomic studies. The study designs are different; nevertheless, there are CpG sites in specific genes repeated in several DNA methylation studies ( Table 7). The T2D biomarkers most studied are CpG sites in the genes ABCG1, TXNIP, SREBF1, and PHOSPHO1. With six [80][81][82][83][84][85], four [80,81,84,86], three [80,84,85], and two [80,83] studies, respectively. The four methylation sites in those genes are in relevant pathways related to T2D. The gene ABCG1 is involved in cholesterol and phospholipid transport and insulin secretion. DNA methylation at ABCG1 has been associated with fasting insulin and HOMA-IR [87,88]. The gene TXNIP is a component of pancreatic β-cell biology, energy metabolism, and cellular redox regulation. TXNIP downregulation protects against obesity-induced diabetes by preventing β-cell apoptosis and preserving β-cell mass [89,90]. The gene SREBF1 is activated by insulin and contributes to dyslipidemia and hepatic steatosis that occurs in obesity and T2D [91]. The gene PHOSPHO1 is considered to be an attractive target for cardiovascular therapy. It has been found that DNA methylation at this gene is correlated positively with HDL levels [83,92]. The results of these studies open the way for functional studies to define their use as biomarkers of T2D and other metabolic comorbidities.
There are other DNA methylation studies of CpG sites that were reported. There were five publications during the 2010-2013 period, the first one discovered three CpG sites in the gene TFAM [93], the second found one CpG in the gene FTO [94], the third found one CpG site in the gene IGFBP-1 [95], the fourth found seven CpG sites in the gene PRKCZ [96], and the fifth found thirteen CpG sites in the gene TCF7L2 [97]. In 2014, an interesting study showed that lower LINE-1 levels after global DNA methylation were associated with a higher risk of a poor carbohydrate metabolic status. This biomarker could be considered a risk factor of T2D and related metabolic disorders, independent of other established risk factors, therefore having a potential role as a possible biomarker [98]. From 2015 to 2017, another four studies suggested new CpG sites of possible interest for T2D biomarker development. The first one suggested five CpG sites in the gene SLC30A8 [99], the second, one CpG site in the gene TXNIP [100], the third one also one CpG in the gene TXNIP [101], and the fourth one two CpG sites, one in the MSI2 gene, and the other one in the CXXC4 gene [102]. At the end of 2017, Van Otterdijk and colleagues reported a tendency between PEG3 DNA methylation levels in peripheral blood leukocytes in T2D patients, and four CpG loci were differentially methylated between groups of participants. As a result, these loci might serve as biomarkers of T2D, although additional research is required to strengthen these observations further [103]. Previous studies had indicated that the ELOVL5 expression is associated with T2D [104,105]. In 2018, Hwang et al. identified the ELOVL5 gene as a novel epigenetic mark in an epigenome-wide analysis of the blood DNA methylation using T2D-discordant monozygotic twins. The results implied that epigenetic alterations in ELOVL5 demonstrate DNA methylation and RNA expression changes. These provide insights into epigenetic biomarkers for T2D [106].
These methylation markers' risk or diagnostic value may be small to justify them as a novel biomarker. Thus, further improvement should include CpG sites showing a higher methylation difference. Most T2D miRNA studies were done on blood samples, predominantly plasma, serum, or exosomes (Table 8). A few studies used whole blood, PBMCs, or endothelial progenitor cells. Almost all studies used qRT-PCR: profiling a large number of miRNAs using PCR panels, targeting a small number of miRNAs, or validating results from microarray, RNA sequencing, or nanostring. Some studies used RNA sequencing without qRT-PCR validation.

DNA Methylation
The discovery of biomarkers in DNA methylation using urine samples in diabetes is currently an interesting topic of research. Cell type-specific DNA methylation patterns have been widely studied and used to calculate proportions of particular cell types related to diabetes in blood but not in urine [173]. Studies focusing on biomarkers for the diagnosis of T2D have not been validated in a clinical environment. Instead, the researchers have been more interested in developing biomarkers for one specific comorbidity, kidney disease (Table 9) [173][174][175]. Because urine examinations identify changes in gene expression in urine-derived cells [176], they may be helpful for the noninvasive assessment of kidney condition and prediction in early T2D [177].

microRNA
T2D urine miRNA studies mainly studied T2D patients with kidney complications [178][179][180][181][182][183][184], except one study compared patients with different T2D treatments [185] (Table 10). Most studies profiled miRNAs from urinary exosomes, while a few profiled urine supernatants [178,184]. There is little overlap between the reported miRNAs because some studies selected miRNAs based on their study population. Nevertheless, most miRNAs were upregulated in patients with T2D complications than patients without complications. The results suggest that urine could be an interesting non-invasive source of T2D biomarkers, but more studies are needed to validate these findings. DNA methylation studies using saliva samples have been more frequent in recent years. Despite this, DNA obtained from saliva is a combination of bacterial and human DNA. Hence, there are unique challenges in using DNA from saliva samples in methylation studies that can affect data quality in epidemiological studies, such as T2D [186]. Currently, there are no empirical studies with evidence for DNA methylation in salivary tissue as potential diagnostic biomarkers of T2D. Nonetheless, some studies support the use of DNA methylation from saliva as a prognostic biomarker in comorbidities. In 2011, 39 genes were identified as differentially methylated in two CpG sites in each gene in a group of end-stage renal disease versus diabetes patients without nephropathy. These sites have been involved in kidney development or diabetic nephropathy. As a result, these sites identified in DNA extracted from saliva give proof that they can be used as predictive biomarkers of disease susceptibility [187]. In 2020 a study provided empirical evidence supporting DNA methylation in salivary tissue as a potential predictor of childhood obesity, making NFR1 a target for further exploration as a biomarker [188].

microRNA
There is only one salivary miRNA study for T2D biomarkers. The study evaluated four miRNAs that have been reported in T2D and periodontal diseases: hsa-mir-146a, hsa-mir-146b, hsa-mir-155, and hsa-mir-203 [189]. These miRNAs were elevated in T2D patients with healthy periodontium compared to healthy controls, with fold changes from 4.9 to 16.5, and were able to distinguish groups with high sensitivity and accuracy (0.75-1 and 0.78-0.85, respectively). The only miRNA that was significantly reduced in T2D patients with periodontitis was hsa-miR-146a, which could predict periodontitis among T2D patients with an accuracy of 0.86. So, salivary miRNAs are promising T2D biomarkers, and more studies are needed to identify more non-invasive biomarkers. Table 11 summarizes and classifies candidate omics biomarkers reported in at least two studies. Most of these biomarkers are risk, diagnostic, prognostic, and response biomarkers profiled from blood samples. To the best of our knowledge, no monitoring DNA or RNA biomarkers have been reported for T2D. So, more longitudinal studies are needed to identify monitoring biomarkers and novel risk biomarkers. Despite many biomarkers reported, there is little overlap between studies due to study methodology differences. For genomic studies, the PRSs can potentially improve biomedical outcomes via precision medicine on T2D. However, large-scale biobanks with diverse ethnic populations need to be considered in more studies to address the need for a broad genomic approach to this global disease.

Conclusions and Future Directions
Biomarker overlap is mainly seen in DNA methylation and miRNA studies (Table 11). Twelve genes were reported in at least two studies for genomics studies and were used as risk biomarkers (Table 11). Only three diagnostic and prognostic biomarkers were reported in at least two mRNA studies (Table 11). The reason is that most mRNA T2D studies are targeted, with researchers focusing on previously reported T2D-related genes. In contrast, most genomics and epigenomics studies used untargeted approaches. Future T2D studies should consider untargeted approaches to identify novel biomarkers.
Lastly, more research on non-invasive biomarkers is needed, especially on T2D patients with no complications. Non-invasive biomarkers would be better for screening and diagnosing T2D in many people. Validation studies are also necessary to confirm the utility of these biomarkers. Nevertheless, these studies have identified potential T2D biomarkers that could help guide future T2D studies.